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Abstract 

We develop algebro-geometrical approach for the open Toda lattice. For a finite 
Jacobi matrix we introduce a singular reducible Riemann surface and associated Baker- 
Akhiezer functions. We provide new explicit solution of inverse spectral problem for a 
finite Jacoby matrix. For the Toda lattice equations we obtain the explicit form of the 
equations of motion, the symplectic structure and Darboux coordinates. We develop 
similar approach for 2D open Toda. Explaining some the machinery we also make 
contact with the periodic case. 
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1 Introduction 



Until now the methods of integration for periodic and open Toda lattice were absolutely 
unrelated to each other. The periodic Toda, is a Hamiltonian system of iV-particles with 
the Hamiltonian 

N 2 N 

H = H^ + Y1 e qk - qk +\ q N+1 = qi + Jo, Io = const. 

k=i z k=i 

The solution of equations of motion was obtained by Krichever using general algebro- 
geometrical approach based on the concept of Baker-Akhiezer function jlj. The solution 
was expressed in terms of theta functions of the spectral curve, associated with the auxil- 
iary spectral problem for an infinite periodic Jacobi matrix. For the open Toda with the 
Hamiltonian 

N 2 JV-1 

H = — + e qk ~ qk+1 
k=i ^ k=1 

the solution of equations of motion was obtained by Moser |J. Moser reduced the original 
problem of integration to the inverse spectral problem for a finite Jacoby matrix. The latter 
was solved by Stieltjes || more than a hundred years ago using continuous fractions. 

Though the Hamiltonian of open Toda lattice can be obtained from the periodic one as 
the limit Iq — > oo, the methods of solution are different. This is due to the fact that auxiliary 
spectral problem associated with a finite Jacoby matrix, unlike the spectral problem for the 
periodic infinite Jacoby matrix, does not determine a natural spectral curve. The first 
attempt to singularize smooth spectral curve of the periodic Toda to obtain a spectral curve 
for the open Toda chain goes back to seventies ||J]. Recently, the interest in this problem 
was revived [|, || due to connections of the Toda lattice with Seiberg-Witten theory of 
supersymmetric SU(N) gauge theory [0. (A relatively complete list of references on new 
insight upon the role of integrable systems in Seiberg-Witten theories |5|-|T^| can be found 

i3i-|3 and books H0-) 
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The spectral curve proposed in [|, |5| for the open Toda lattice is determined by the 
equation 

N 

w Q = P(E) = H(E-E k ) (1.1) 

k=l 

and can be considered as the limit when I Q — > oo of the hyper-elliptic spectral curve 

w + — = P(E), A 2 = e- /o (1.2) 

w 

for the periodic case. 

From algebro-geometrical viewpoint the limit A — > of hyper-elliptic curve leads to a 
singular curve, which is two copies of rational curve, corresponding to two sheets of hyper- 
elliptic curve, glued together at points E^. It is well known that Baker-Akhiezer functions 
introduced originally for smooth algebraic curves, are also an important tool of integration 
in the limiting case of singular curves. Multi-soliton and rational solutions of integrable 



2 



equations can be obtained within this approach (see in [19|). The main goal of the present 



paper is to demonstrate that algebro-geometrical approach based on the concept of Baker- 
Akhiezer function can be used in the case of reducible singular algebraic curves. As an 
outcome, we provide a solution to the inverse spectral problem for a finite Jacoby matrix 
which is different from the classical Stieltjes' solution. 

In the recent papers by Krichever and Phong, |2(| [21], the new approach using Baker- 
Akhiezer functions for construction of Hamiltonian theory of soliton type equations was 
developed. It provides a universal scheme for construction of angle-action type variables for 
these equations. In Section 4 we illustrate these ideas, and show that, the case of the open 
Toda lattice with singular curve can be treated similarly. 

Finally, in section 5 we indicate how the ideas introduced in this paper can be applied 
to integration of open two-dimensional Toda lattice. 



2 Periodic Toda lattice. 

To begin with, let us first briefly recall algebro-geometric solution of the periodic Toda lattice. 
That will help us later to clarify algebro-geometric origin of our approach to the open Toda 
lattice. 



Most of the material is standard and details can be found in 22, uM- It is convenient 



to consider the periodic Toda lattice system as a subsystem of an infinite lattice. As it was 



found in p| , pq] , the equations of motion 

q k =Pk, Pk = -e« +1 +e qk -'- qk . (2.1) 

k = —1, 0, 1, ... are equivalent to the Lax equation L = [A, L] for auxiliary linear operators, 
where L and A are difference operators 

(Lif)) n = c n ip n+1 + v n ip n + c n _iVv-i (2.2) 

{M)n = yi-n+l - -y^n-l ( 2 -3) 

with coefficients 

Ck = e (« + i)/2 ; Vk = _ pk _ (2.4) 

Let q n satisfy the constraint q n + h = Qn+N which is invariant with respect to (^j]). The 
corresponding L and A operators are periodic ones: c n = c n+ N, v n = v n+ N- The Floquet- 
Bloch solution is the solution of the periodic Schrodinger equation 

(Llp) n = C n 1p n+1 + V n 1p n + C„_i^„_i = Elp n (2.5) 

such that 

ip n+N = wip n . (2.6) 

Let T(E) be a restriction of the monodromy operator {Tf) n = f n+ N on the invariant space of 
solutions to the equation (Lip) n = Etj) n . The solutions if and 6 of the Schrodinger equation, 
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normalized such that (po = 1, y?i = and 6 = 0, Q\ = define a basis in which the 
monodromy operator have the form: 



T(E) 



9 N 
9n+i 



(2.7) 



Then pairs of complex numbers (w, E) for which there is a common solution of equations 
( p.5| , p.6|) are defined by characteristic equation 



R(w,E) = det(w-T(E)) 



w 



2A{E)w + 1 = 0. 



{21 



The hyperelliptic Riemann surface Y defined by equation fl2.8|) is called a spectral curve of 
the periodic Schrodinger operator. A point on the curve we denote by Q = (w, E). 

To make contact with the open Toda case it is useful to introduce a different repre- 
sentation of T. Consider the operator L(w) which is defined as a restriction of the infinite- 
dimensional operator L on iV-dimensional space of functions that satisfy equation (|2.6|) . The 
corresponding matrix has the form 



L(w) 



From (EOT) it follows that 



Co 




WC N _i 



Co 





Ci 

CjV-3 v N-2 

Cat„ 2 






CN-2 

v N -i 



(2.9) 



det(E - L(w)) = P(E) -Aw- Aw' 1 = 0, A = (-1 
Coefficients of the monic polynomial P{E) 

P{E) = E N 



N 

Ii C n 

n=l 



: _l)iV e -7 /2_ 



(2.10) 



N-1 
i=0 



(2.11) 



are polynomial functions of the coefficients c n , v n of L. For example, 



un-i 



N 

E 

n=l 



N 



N 



V n , U N - 2 



E ^-E c i 



(2.12) 



v0<i<j 



i=l 



Note, that P = 2AA(E). Change of the variable = Aw transforms ( |2.10| ) into the desired 
formula (O). 



The equation d w R(w, E) = for branch points of T leads to A 2 (E) = 1. It has 2N 
real roots which are points of periodic/anti-periodic spectrum for L, according to the sign of 
A(E) = ±1 at these points. For a generic configuration of particles these roots are distinct 
Ei < E 2 < ■ ■ ■ < E 2 n and the curve V is smooth and has genus N — 1. 

To make a model of the curve V we take two copies of the complex plane and make N 
cuts along the bands [Ei, E 2 ], ■ ■ ■ , [E 2 n-i, E 2 n]- Then glue two cut-planes together. 
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The multiplier is holomorphic in the affine part of the curve with pole/zero of degree N 
at the infinity P + /P_. The multivalued function of quasimomentum p(Q) is introduced by 
the formula w = e Np . It has the asymptotic expansion 

OO TT 

±p(E) = log E E = E(Q), Q-(P±). 
The coefficients are standard integrals of the periodic Toda 

T 1 N 1 N „2 N 



The explicit formula for the Floquet-Bloch solution 



k- 



MQ) = ME) + W{Q) -^ {E) e n (E). (2.13) 

t>N{E) 



implies that 



Lemma 1 (i) The Floquet-Bloch solution, normalized by the condition ip — 1 becomes 
single-valued on the Riemann surface T. It has a single pole jk on each real oval a^, k = 
1, . . . , N — 1, which is a preimage of the [E 2 k, E 2 k+i\- 

(ii) In the vicinity of infinities the Floquet solution has the asymptotics 
ME) = E ±n e^-^l 2 (l + pxt(n)E-'j , 
where E = E(Q), Q -> (P±). 

Equation ( |2.13| ) implies that projections E{^ s ) of the poles of ip are zeros of the polynomial 
9n(E) : 9n(E( , j s )) = 0. Therefore, 9 n (E(j a )) satisfies zero boundary conditions, and E{^ s ) 
are points of the Dirichlet spectrum for the Schrodinger operator. 

The map which associates for iV-periodic operator L with coefficients {c n ,v n } a spectral 
curve T, and the divisor D = {ji, . . . ,7at_i} of poles of the Floquet-Bloch solution. 

{c n ,v n } -> {r,P}, (2.14) 

is referred as direct spectral transform. 

The spectral curve is defined by the variable A , and iV coefficients Ui of the polynomial 
P{E). Therefore, a space of the spectral data has the same dimension 2N as the space of 
iV-periodic operators. It is a fundamental fact in the theory that the map ([2.14|) is a bijective 
correspondence of generic points. We shall refer to the reverse construction 

{T,D}^{c n ,v n }, (2.15) 

which recaptures the dynamical variables {c n ,v n } from the geometric data {T,D}, as the 
inverse problem. As usual in the algebro-geometric theory of solitons, it will be based on the 
construction of a Baker- Akhiezer function. 
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Lemma 2 /fj/ Let V be a smooth hyperelliptic curve defined by equation (\2.1Q) . Then for a 
generic set of (N — 1) points 7 S G T there exists a unique (up to a sign) function ^ n (t,Q) 
such that: 

(i) it is meromorphic function on V outside infinities P± and has at most simple poles at 
the points 7 S . 

(ii) In the vicinity of infinities the Floquet solution has the asymptotics 

9 n {E) = E ±n e ±Et '%f{t) (l + jTxf(n,t)E~ s ) , 
where E = E(Q), Q -+ (P±). 



Uniqueness of \1/ immediately implies the following theorem. 

Theorem 1 The Baker- Akhiezer function ^ n (t,Q) associated with the spectral data {T,D} 
satisfies equations 

(LV) n (t,Q) = E* n (t,Q), (dt-A)* n (t,Q) = 0. (2.16) 



where the operators L and A have the form (2.i ) and ( \2.3j) with coefficients defined by the 
formulae 

c„(t) = v n = xt(n,t) - xt(n + l,t). 

These coefficients are N -periodic functions of the variable n. 

From the Lax equation, which is a compatibility condition for ( |2.16| ) we obtain that 

q n {t) = \nb n {t) (2.17) 

is a solution of the periodic Toda lattice normalized by the condition go(0) = 0. Note, that for 
t = the Baker- Akhiezer function has the same analytical properties as the Floquet-Bloch 
solutions. Therefore, ^/ n (0,Q) = ip n (Q), and using the spectral data {T,D} associated with 
{c n ,Vn} w e solve the Cauchy problem for the periodic Toda lattice. 

This solution can be written explicitly in terms of the Riemann ^-function associated 



with matrix of 6-periods of a normalized holomorphic differentials [22] . 



. x n , fe((n- l)U + Vt + Z)6(Z)\ 
«.(«) = *(0) + jfh + in ( ^ + > vt + z)e[z >_ { v > ) (2.18) 

Here U and V are (N — l)-dimensional vectors that are periods of certain meromorphic 
differentials on T, and Z is Abel transform of the divisor D. 

We finish this brief account of solution of the periodic Toda lattice by the following 
remark. The Lax equation implies that if the dynamical variables {c n (t),v n (t)} evolve ac- 
cording to the Toda lattice equations then the spectral curve of the Schrodinger operator 
is time-independent. At the same time the poles, of the Floquet-Bloch solution normalized 
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by the condition i[) = 1 do depend on t. From fl2.16|) it follows that ipnit^Q) and the 
Baker- Akhiezer function are proportional to each other. Therefore, 

ip n (t,Q) = ^ n (t,Q)%\t,Q) 

Hence, the divisor D(t) can be identified with the divisor of zeros of the function ^o(t,Q). 
An image of this divisor under the Abel transform evolves linearly, Z(t) = Z + Vt. 



3 Open Toda lattice 



Equations of motion for the open N-particle Toda lattice 



Qk 
Pk 
Po 



Pk, 

g<?fc— <?fc+l _|_ g9fc-l-9fc 

oio-qi 



k = 0,...,N- 1, 



l,...,N-2, 



— e* 



Pn-i 



D qN-2—qN-i 



have Lax representation d t L n = [A n ,L ] where L and A are finite-dimensional Jacobi 
matrices 



c 

c Vi 

••• 
••• 





ci 








CiV-3 Vn-2 Cn-2 
C N _ 2 VN-1 



An 



c 

-c 

••• 

••• 





ci 

-cn-s 







C N -2 
-C7V-2 



(3.1) 



where the variables c n , v n are defined by formulae ( |2.4j ), as in the periodic case. 

Equations of motion are invariant with respect to the shift v n — > v n + const. Therefore, 
without loss of generality we will assume from now on that v n are normalized by the constraint 



N-l 
n=0 







(3.2) 



As it was mentioned in the introduction, the spectral problem for L does not contain 
spectral parameter. Let us introduced another auxiliary spectral problem, which does contain 
spectral parameter. The very same equations of motion are equivalent to the Lax equation 
for the following finite-dimensional operators: 



L(w) 



v 
Co 





Co 

Vl 





Cl 
C7V-3 








v N-2 CN-2 
Cn-2 UJV-I 



(3.3) 
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A(w) 



c 
-c 





Cl 








where 





WC N _i 

Cjv-i = e 



-C N - 3 Cat_2 

- Cjv _ 2 



(3.4) 



(giv-i-go)/ 2 



AT-2 



n c- 



(3.5) 



n=0 



Coordinates ip n (E), n = 0, . . . , iV — 2, of an eigenvector t/j, L(w)ip(E) = Eip(E) can be found 
reccurently from equations 



Eip = coV'i+^oV'o, 

Eip n = c n ip n+1 + v n ip n + c„_i^„_i, n = 1, . . . , N - 2. 
If we take ?/>o = 1, then fl3.6|) implies that i/j n (E) is a polynomial of degree n 



iP n (E)=J2Hn)E j , b n (n)=e 



(9n-9o)/2 



i=0 



The last equation 



£?ViV-l = C N -iW1pQ + WAr_i^Ar_i + C N - 2 ^N-2 



(3.6) 



(3.7) 



(3i 



defines w = w(E). It is a polynomial of degree iV. From ( |3.5| ) it follows that w(E) is a 
monic polynomial. It can be also found from the characteristic equation 



N-2 



Rq(w,E) = det(E-L(w)) = P{E) - w = 0, P{E) = E N + £ UiE\ (3.9) 



i=0 



which determines the components of the spectral curve. For w = the matrix L(w) coincides 
with Lq. Therefore, zeros Ek of the polynomial w(E) 



N 

w{E) = 1[{E - E k ) 

k=l 



(3.10) 



are eigenvalues of the matrix Lq. The spectral curve To is singular algebraic curve which is 
obtained by gluing to each other two copies of the E'-plane along the set of N points Ek- 

Let us introduce now a solution (p = {(p n (E)} of the adjoint equation L T (w)(p = Eip. If 
we normalize <p by the condition (p^-i = 1, then the equations 



Eip n = c n ip n+1 + v n ip n + c n _iv?„-i, n = N-2, 



1. 



reccurently define (p n (E), n = 0, . . . , iV — 2. It is a polynomial of degree iV — 1 — n with the 
leading coefficient 

<p n (E) = e (lN-l-qn)/2 E N-l-n + Q^N-2-n^ 
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The last equation 

Eif = c (fi + two + wc N -i<p N -x, 

which determines w leads to the same formula ( [3.10|) . 

Let j s be roots of the equation (po(E) = 0. As we will show, they interlace N roots of 
the determinant of Lq, i.e. 



E x < 71 < E 2 < . . . < £tv-i < 7jv-i < E 



N- 



(3.11) 



Coordinates of an eigenvector ip a (E) = ipn(E) for the matrix L T (w), normalized by the 
condition ip% = 1, are equal to ip* = </V 1 (E)(p n (E) , and therefore, have the form 



E? = -Q- n b°(n)Ei 



3 (90-9n)/2 



(3.12) 



At w = matrices L(w) and L T (w) coincides. Therefore, their eigenvectors are proportional 
to each other. Due to our choice of normalization, we conclude that 



ip n {E k ) = r n {E k ). 



(3.13) 



Due to (|3.13|) , a pair of functions i/j n (E), if)£(E) can be considered as a single- valued function 
on To- This is the Baker- Akhiezer function for the spectral curve IV 

To prove ( |3.11|) we represent the solution <p using truncated matrices 



v k c k 













CN-3 VN-2 Cn-2 
C N _2 V N -i 



(3.14) 



We define det{E - L N+1 ) = 0, 
det(E - L k ) = 
it is easy to check, that 



det(E — L N ) = 1. Using the reccurence relations 
(E - v k ) det(£ - L k+1 ) - c 2 k det(E - L k+2 ), 



det(E - L n+1 ] 



C n ....CN-2 

From this we see that 7 S are the spectrum of the matrix L\ which is obtained from L by 
canceling the first raw and the first column. Now classical theorem of Sturm |26j implies 
( p. 11] ). By analogy with the periodic case we call the points j s the divisor D . The points of 
D Q move between E k under the Toda flow, but never change a sheet, contrary to the periodic 
case. 

As in the periodic case, a map which associates to dynamical variables a set of spectral 
data Tq, -Do of the same dimension is a bijective correspondence 



{C , . . . , Cjv_ 2 , V , ■ ■ ■ , WjV-l} < — > {IV D }. 



(3.15) 
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In fact, the gluing conditions ( p. 13 ), and the normalization 

= 1, (3-16) 



which follows from (|3.7|) and ( |3.12|) , allow to reconstruct ip n and if>°- 



For each n the gluing conditions are equivalent to N linear equations for N + 1 unknown 
coefficients bi, b". If we define the vector with N + l components Bj = bj, j — 0, . . . ,n; Bj = 
bj-n-v 3 — n + 1, . . . , N then the matrix of the corresponding linear system has the form 

Mij(n) = nEl j = 0,...,n 

M l3 (n) = -Et n '\ j = n + l,...,N, (3.17) 

where = H^x(Ei — 7 S ), i = 1, . . . , N. Let Aj(n) be the j—th. minor obtained by canceling 
the j-th column. The coefficients B's which satisfy the normalization condition are 

/[,{„) , ■ ./ V. (3.18) 



JA n (n) A N (n 

Since B n (n) = e (9 ' l - 9o)/2 , then 



A (n) 

q n -q = 21nB n (n)=ln-^r, n = 0, . . . , N - 1. (3.19) 

L\ N {n) 



Of course, there are sets of complex numbers {Ek, 7 S }, such that the right hand side of ( |3.19| ) 
is singular. 



Lemma 3 Let the spectral data {Ek, j s } be real and satisfy conditions ( 3.11 ), then 



A n (n)Ajv(n) ^ 0. (3.20) 

Indeed, let ip n and be the unnormalized solution of the gluing conditions with the coef- 
ficients bi, bj, corresponding to Bj(n) = Aj(n). The product ip n ip^ * s a rational function of 
the form 



W» = TW^7§ ] ^ » F = An(n)A^(n)^- 1 + 0(E N -^ 
1L=i (E -j s 



If the spectral data are real, then ip n (E) and if}^{E) are real for real values of E. Due to 
the gluing conditions the product of these functions at E = Ek is positive. Therefore, the 
product ipnipn has at least one zero on each interval [Ek, Ek+i], because it has a pole there. 
Hence, degree of the polynomial F can not be less then N — 1. That implies ( |3.20| ). 

Theorem 2 Let ip n {E),ip^{E) be functions associated with spectral data Ek,j s . Then they 
satisfy equations L(w)ip = Eip and L T (w)ip a = Eip a , where 

, A n (n)A N (n + 1) 
A n+1 (ri + l)Ajv(n) 
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and 



Vn = " u / - " v -I '— n=l,...JV-2; 
A n (n) A n+1 (n + l)' 

A (l) Ajv- 2 (iV-l) 



A!(l)' ^ A^iV-l)' 
Proof. Functions ip n , ifi* at E — > oo have the form 

Vn(£) = EV^Wl + f] V (3.22) 



s=l 



<(£) = E' n e {q °- qn)/2 h+^x«W^" s J • (3.23) 

Note, that 

*.(») = XTW = ^4| + E7,. (3.24) 

Here for n = and n = N — 1 we formally put 6-i(0) = b a _x{N — 1) = 0. 
Let 4> n (E) and 0£(-E) be coordinates of the vectors 

0(E) = (LH - E)^(E), r(E) = (L T (w) - E)r(E). (3.25) 

If c n are defined by (|2.4|), and u n are defined by the equations 

«n = Xi(n)-Xi(n + l), n = 0,...,iV-2; = X i{N - 1), (3.26) 

then ( |3.22| ) implies that near infinity has the form <p n = 0(E n ~ 1 ). Therefore, we conclude 
that (p n is a polynomial of degree at most (n — 1). The leading term of expansion for 0^ near 
infinity is 0(E~ n ). Hence, 0^ has the same form as Functions n and 0^ satisfy the 
same gluing conditions ( |3.13| ). They are equivalent to N linear equations for N unknown 
coefficients for corresponding polynomials. The matrix of this system is obtained from the 
matrix M ( 3.171) by canceling the n-th column. If this matrix is non degenerate, A n (n) ^ 0, 



we conclude that = 0°" = 0. Formulae ( |3.18| , |3~2"4|) allow to complete proof of the theorem. 
Note, that the equation (L T — E)ip a = implies that 

vo = -xl(N-l), v n = x a x(n)-X a x(n-l), n = l,...,N-l, (3.27) 
Therefore, we obtain the identities 

Xi{n + l)+xi{n) = -v , (3.28) 
which will be used in the next section. 
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Solution of the Cauchy problem for the open Toda lattice can be obtained in a way 
which is almost identical to the periodic case. First of all, we defined time-dependent Baker- 
Akhiezer functions ip n (t, E), if)%(t, E), as functions of the form 

Mt,E) = e Et / 2 (^h{t,n)EiY (3.29) 

= ^ ^g ), ,3,0) 

that satisfy gluing Q3.13|) and normalization conditions (|3.16|) . 

After that, without any change of arguments we obtain that they satisfy the equations 
{Lit, w) - E)ip(t, E) = 0, {L T (t, w) - E)tp a (t, E) = 0. In the similar way we get that these 
functions satisfy equations (d t — A(t,w))ip(t,E) — 0, (dt — A T (t,w))ip a (t,E) =0. As a result 
of compatibility condition we get that functions q n (t) = qo{0) + lnb n (t,n) are solutions of 
the open Toda lattice. 

The matrix of the linear system for bi(t,n) and 6j(t,n) has the same form as in ( |3.17| ) 
with the only difference that constants r; should be replaced by rj(t) = rie Ekt . Therefore, 
we come to the following theorem. 

Theorem 3 Let Aj(t,n) be the j-th minor obtained by canceling the j-th column of the 
N x (N + 1) matrix 



M i3 (t } n) = riE{e Ei \ j = 0, 



n 



Myfan) = -El~ n -\ j = n+l,...,N, 

where = Y\^=i(Ei — 7 S ). Then the formula 

A n (t,n) 



q n {t) = 9o (0) +ln- 



A N (t,n) 

gives a solution of the Cauchy problem for the open Toda lattice equations. 



4 Action-angle type variables 

The main goal of this section is to show that Lax representation which depends on the 
spectral parameter allows to apply for Hamiltonian theory of open Toda lattice an algebro- 
geometric approach proposed in |2D], |2"I| . The main idea of this approach is to introduce in 
a universal way two-form on a space of auxiliary linear operators, written in terms of the 
operator itself and its eigenf unctions. This form defines on proper subspaces a symplectic 
structure with respect to which Lax equation is Hamiltonian. Moreover, the way to define 
the symplectic structure leads directly to construction of Darboux coordinates. 

In the case of consideration, we define such a form by the formula 



L) = ies r 



d In xv 

< xp + SL A dip > + < 54> + A 5L4> > . (4.L 
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Here 5L is an operator-valued one-form on a space of coefficients of the operator, i.e. in our 
case it is just a matrix with entries, which are one-forms 



SL(w) 



5v 5c q 
5c 5v\ 5ci 





u>5cjv-i 










5v N _2 5c N _2 
5c N „2 <^jV-l 



(4.2) 



An eigen-vector ip of L can be seen as a vector-valued function on the space of the coefficients 
of L, and therefore, its external differential Sip is a vector-valued one-form. The co(row)- 
vector ip + is an eigen-vector of L in the dual space, ip + L = Eip + , i.e. ip + = {ip cr ) T - At last, 
< f + g > stands for pairing of vector and co- vector 



N-l 



< f + 9 >= E fn9n- 



n=0 



Now, it is necessary to add an important remark, which clarifies our understanding of 
dip. In ( |4.2[ ) variation of L was defined such that the spectral parameter w is fixed. An 
eigenvalue of L can be seen locally as a function E = E{w) defined implicitly by ( |3.9|) . 
Therefore, though ip(E) is a polynomial in E, locally, we consider it as a function of w, 
ip(w) = ip(E(w)), and partial derivatives in 



8 



E fa 



_d_ 

OCr, 



SVr, 



d_ 



are taken for fixed values of w. A connection between 5 and variation <5 , which is defined if 
one keeps E fixed, is just a chain rule. Let f(E) = f(E, c n , v n ) be a function of the variable 
E depending on {c n ,v n } as on external parameters. Then 



5f(w) = 5f(E(w)) = 5 f(E) + 4L8E(w). 



(4.3) 



In particular, if we take f = w, then from (O) it follows that 



= 6 Q w(E) + ^8E(w). 



Therefore, ( f4.3| ) may be rewritten in the form 



df 



5f(w) = 5 f(E) - ^ow(E). 



(4.4) 



This formula shows that variation 5 leads to an appearance of poles at zeros of dw , which 
are critical points of the polynomial w(E), where we, even locally, can not invert w and 
introduce E = E(w). 



Our immediate goal is to find an explicit formula for uj in term of dynamical variables. 
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Lemma 4 Let to be given (J71 ). Then on symplectic leaves of the space of dynamical vari- 
ables {c n ,v n }, defined by the constraint J2n v n = const, the form u is equal to 



N-l 



it) 



2 ^ <5g„ A 5v n 



(4.5) 



n=0 



Proof. In order to simplify slightly formulae, we consider only the leaf J2n v n = 0. For all 
the other leaves a proof is going in the same way. 

At infinity the functions ip n and ip a can be expanded as Laurent series in the variable 



From (|S^J5^ ) we get 



ip n (k) 



k = w 1/N {E) = E + 0{E- 1 ). 

(oo 

/ oo 

= k- n e {q °- qn)/2 1 + E £(n)k- 



(4.6) 
(4.7) 



=i 



N -2, 



(4.8) 
(4.9) 



Note, that in this expansion £i = xi, £i — Xi- Therefore, as before, 

v n = 6 (n) - £i(n + 1), + 1) + = -v n = 0, 

£ l (N-l)=v N - 1 . 

The variable A; is a constant with respect to our definition of 5, and therefore, an expansion 
of Sip can be obtained by taking variation of coefficients in ( |4.6| , f4.7|) . 

Definition of 5L implies 

N-l N-2 



< 



n=0 



n=0 



From (4.S , 4.9) it follows that at infinity 

din w 



Using expansions (|4.6j , f4.7|) and equation (|4.8|) we obtain 



(4.10) 



u = -2 



7V-1 JV-2 

E 6v n A 5(g n - go) + E ~~ Qn+i) A 5£i(n + 1) 

n=0 n=0 



(4.11) 



due to the relation ^ n = an d equations ( |4.8|) . The lemma 



This equation implies 
is proved. 

Representation of symplectic forms as a residue at infinity of the spectral curve of mero- 
morphic differential dQ in the right hand side of ( |4.1[ ) allows to find Darboux coordinates. 
These Darboux coordinates are connected with residues of the differential at finite points. 
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Lemma 5 Let R = Yl s (E — 7 S ) be the monic polynomial with zeros at the poles of the 
Baker- Akhiezer function ip a . Then the following identity holds 



dw 



< ip+ip > 



RdE, 



(4.12) 



Proof. The equation (L — E)ip = implies that (dL — dE)ip = (E — L)dip. Therefore, 

< ip + (dL - dE)ip >=< ip + (E - L)dip >= 0, 
and we obtain the equality 

dE < ip + ip >=< ijj + dLijj >= CN-idw (Vtv-iV'o) • 

From (|3.12| ) we have that V'jv-i = C/v-i-^ -1 • Lemma is proved. 

This result implies that the differential dQ, besides E = oo, may have poles at spectral 
points Ek, at poles 7 S of ip + , and at zeros Zi of dw. (The last poles are due to ( |4.4|) .) 

First of all, let us show that dfl has no poles at E^. At these points L = L , and therefore, 
is symmetric. Using gluing conditions we obtain 



< i) + 5L A 5ip > + < 5i/j + A 5Lijj > 



w=0 



< ip 1 SL A Sip > + < Sip 1 A SL ip >= 



due to skew-symmetry of the wedge product. 

At E = 7 S the first term of dfl does not contribute to residue because poles of ip + cancel 
with zero of R. Nontrivial residue is due to second order pole of Sip + . We have 



res^ 



< Sip + A SLip > 
< ip+ip > 



dlnw 



<51nu>(7 s ) A 



< ip + SLip > 



E=-y s 



(4.13) 



< ip + ip > 

= S\nw(^s) ASE^s). 

In the last equality we use the well-known formula for variation of eigenvalue 

< ip + SLip >= SE < ip + ip >, 

which immediately follows from the equation 

< ip + (SL - SE)iP >=< iP + (L - E)Sip >= 0. 

Consider now critical points Zi of the polynomial w(E) : dw{zP) = 0. Formula 
implies 

<ip + SLASip> „ [< ip + SLdip > ASE(w)} dlnw 

res 2 . — — dlnw = res z . — — — 

1 < ip+ip > dE < ip+ip > 

Due to skew-symmetry of the wedge product, SL in the last formula may be replaced by 
(5L - SE). The equations ip + {SL - SE) = -Sip + (L - E) and (L - E)dip = -(dL - dE)ip 
imply 

<iP + (SL- SE)diP >=< Sip + (dL- dE)iP > . 



for Sip 



(4.14) 
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Note that dL(zj) = 0. Hence, 

< ip + 5L A 5ip > „ (< 5ib + ib > A5E) dlnw 
res 2 . — — dlnw — —res,- — — . 

1 < t/)+ip > < ip+ip > 

In the same way we obtain 

< 5ip + A 5Lijj > „ (5EA < ip + 5ip >) dlnw 
res z . — — ■ amw = — res z . — — . 

< ip+ip > < ip + ip > 

Taking a sum of these equalities we get 

(< ip + 5ip - 5i/j + iIj > A5E) dlnw 

res,. cm = res,. ; . 

< ip+ip > 

Besides poles at the points z^, the differential at the right hand side of the last equation has 
poles at the points 7 S , only. Indeed, it has no pole at the infinity, because due to constraint 
52 n 8v n = we have 5E = 0(k~ 1 ). There is no pole of the differential at zeros E k of w(E) 
due to gluing conditions. 

Hence, we get 

£ res,,* = - £ res. « t*t - > -gS "■"» = £ 6 In u>( 7> ) A SE M (4.15) 

Sum of all the residues of dfl equals to zero. Therefore, using (|4.13| , |4.15|) we obtain 
Lemma 6 Let uj be a two- form given by the formula (U-A) - Then 

N-l 

u = 2 SE( ls ) A 5 In w{n s ). (4.16) 

8=1 

Finally, combining results of Lemma 4 and 6 we obtain 

Theorem 4 Poles 7 S = E(n s ) of the Baker- Akhiezer function and evaluation lnu>(7 s ) are 
Darboux coordinates for reduction of the canonical symplectic form on the leaf Mi, defined 
by the constraint J2n v n — d: 

\y^5v n A5q n \ = E Slnw(j s ) A 5^ s . 

\ n J Mi s=1 

From (|3.9| ) we have iw(7 s ) = llfc(7s — Ek)- Therefore, 

N-l N-l N srp a s N 

E 6hxw{>y.) A 5 7s = E E Z" _ = E^ ln ^ A SE k . (4.17) 

s=l s=l k=l ^ k ^ s k=l 

where r k = R(E k ) = Yl s (E k — j s ) are variables used in the previous section. From the 
formulae ipo{E k )ip N -i{E k ) = 1 and <p (E k ) = CN-\r k) k = 1, ...,N; we obtain 

N N N 

J2§lnr k ASE k = S ^ ip {E k ) A 5 E k = - E $ ln^-i(-Efc) A SE k 

k=l k=l k=l 
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Note that in all formulae the variables are subject to the constraint J2k^k = J2n v n — I- 
The last representation for the symplectic structure in terms of first/last components of 
normalized eigenfunctions and eigenvalues of L is well-known, see [f^]. 

We would like to emphasize that the form in ( 4.171) can be seen as a limit A —>■ of the 
formula for action-angle variables for the periodic Toda lattice. To see this, let us choose 
Ei, . . . E N -i, as independent variables, then 



N N-l 

J2^\nr k A5E k = ^ 5(\nr k -\nr N ) A 5E k . 

k=l k=l 



(4.18) 



If we choose first (N — 1) cuts of the hyperelliptic curve as a basis of a-cycles, then E k are 
identified with limit of the action variables for the periodic Toda lattice: 



E k = lim / EdlnnjQ 



where wo = Aw from (|1.2|) . Limit of the normalized holomorphic differential dVt k on the 
spectral curve of the periodic Toda lattice is equal to 



lim dQ k 



dE 



dE 



(E - E k ) (E - E 



Therefore, the variables (lnr^ — lnrA?) are identified with a limit of the angle variables <p k 
in the periodic case, which are result of the Abel transform of poles of the Floquet-Bloch 
solution 



In r k — In tn 



lim 



lim 

A— 



5 2D open Toda lattice 

Equations of motion of 2D open Toda lattice are equations for unknown functions q n (x, t). 
In light-cone coordinates £ = x + 1, rj — x — t they have the form 



gQn—qn + 1 



1,...,N 



oQN-2—QN-l 



(5.1; 



We consider the Goursat initial value problem with initial data given on characteristics 
%n(£, 0) and q n {0,7)). 

In two-dimensional case an analog of the Lax representation is the zero-curvature equation 
— U , d v — Vq] =0, where U and V are finite-dimensional matrices 



U 



v 1 
v 1 1 








•■•0 v N - 2 1 
•■• v N -x 



V 





cl 

••■ 

••■ 



c 






2 

N-3 





... o 
... o 





(5.2) 
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where, as before c\ = e Qn qn+1 and v n = —d^q n . 



Note, that for N = 2 equations (|5.1j ) imply for <p = q — q x the Liouville equation 

%fp = -2e*. 

Explicit solution for the Liouville equation which contains two arbitrary was constructed by 
Liouville himself ||27f . There are many approaches that lead to generalization of this solution 
for the case of arbitrary N (see, for example, |28|. |29| ). 

The main goal of this section is to show that the general solution of the open 2D Toda 
lattice can be obtained as a degenerate case of the construction proposed in [|30| for solving 
of periodic 2D Toda lattice. In j30[ it was shown that any local solution of periodic 2D 
Toda lattice has an analog of the d'Alambert representation as a composition of functions 
depending on £ (or rj), only. Nonlinear superposition of such right- (or left-) moving waves 
is achieved with the help of auxiliary linear Riemann-Hilbert problem. The corresponding 
matrix Baker-Akhiezer function which solves the Riemann-Hilbert problem has essential 
singularities at two points on complex plane of the spectral parameter w. If we represent 
this plane as two-fold cover with the help of the equation w + A 2 w~ 1 = E and take A — > 0, 
then the limit of the two-fold cover becomes a singular curve, which can be seen as two copies 
of the plane attached to each other at one point. In this limit, as we show the Riemann- 
Hilbert problem reduces to a system of linear equations and provides the general solution 
for 2D open Toda lattice. 

Now let us provide details. Initial data which define the general solution are arbitrary 
functions of one variable ao(£)> ■ ■ ■ > a N-i(Q, bo(r]), • • • , ^-2(^7) which can be identified with 
initial data for q n on the characteristics a n (£) = t> n (£,0), b n {rf) = c^(0,rj). 

Let $(£,«;) and § u (rj,w) matrix solutions of ordinary differential equations 



(3f - w)) w) = 0, (d v - B(<n, w)) $'(77, w) = 0. 
normalized by the initial condition 



(5.3) 



$(0,w) = $ ff (0,u;) 



(5.4) 



Here A and B are the matrices 



a 1 
ax 1 





w 



d]\j—2 








1 

ajv-i 



B(r],w) 




b 







■■• &jv_3 

••• b N - 2 



(5.5) 



As before, we define bN-i by the formula 6jv-i = II„=o u n 

Let R(£,rj) and R a (^,ri) be lower- and upper- triangular matrices 



Ri. 



1, Rij = 0,i<j; 



0, i > j- 
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They are uniquely defined by the gluing condition at the point w = 

77,0) = 77,0) (5.6) 

for matrices 

(£, 77, w) = R(£, 7i)$(£, w), 77, w) = , 77)^(77, «,), (5.7) 

If diagonal matrix H = Hidij, upper- and lower- diagonal matrices B + , B- with units on 
diagonals, are defined by the Borel decomposition of the matrix 

<f>(Z,0)$' T (r],0)- 1 = B_HB + , (5.8) 

then (|5.6|) implies that R = Bz 1 , R a = HB + . In particular we have R^ = Hj. 

Theorem 5 The functions 

q n (^r ] ) = -\nR^v) 
solve the open 2D Toda lattice equations. They satisfy the initial conditions 

%n(£,0) = -a n (C), gn(0,77) - g n+ i(0,?7) = lnb n (rj). 

Proof. Let us show first, that the matrix functions \I/ and , satisfy linear equations 

- U)V = 0, - U )*° = 0, (5.9) 

where Uq = Uq(^,tj,w) has the form ( |5.2|) with v n = —d^q n and U = Uq + we~ . Here e~ is 
the matrix with the only non- vanishing entry at the left lower corner: e~- = S^n-iSqj. 

The matrix \1/ is non-degenerate. Therefore, its logarithmic derivative U = d^^fr" 1 is 
holomorphic function of w. Moreover, from (|5.7|) and condition that R is a lower-triangle 
matrix it follows that 

U = R^R' 1 + RAR- 1 (5.10) 

have the form U = we" + U, where U does not depend on w and have vanishing entries over 
the first diagonal over the main diagonal 

U hl+1 = 1, U lj= 0, j > i + (5.11) 
On the other hand, from (|5.7|) we have that 

U = = R a (R^y 1 (5.12) 

is upper-triangular matrix with diagonal elements equal —d^q n . Gluing conditions ( |5.6| ) imply 

U(&ri,0) = U(&Ti) = U o (£,Ti) (5.13) 
Therefore, combining the conditions (|5.11 ) and the fact that Uq is upper-triangular we con- 



clude that Uq has the form ( pT2 
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In the similar way one proves that ^ and satisfy linear equations 

(d v - V )V = 0, (9 e - V)^ = 0, (5.14) 

where V = Vo + u>(e~) T . Compatibility conditions of ( |5.9| ) and ( |5.14j ) imply that q n solves 
2D Toda lattice equations. 

At rj = due to ( |5.4| ) we have that matrices i?(£, 0), .R s (£, 0) are defined by the Borel 
decomposition of the matrix 

$(e,0)= J R- 1 (e,0) J R ff (e,0) 
But this matrix is upper-triangular. Therefore, i?(£, 0) = 1 and 

tf o (£,0) = A(£). 

In the same way we prove that Vo(0, 77) = B(rj). The last two equalities gives initial conditions 
on the characteristics for q n . The theorem is proved. 

We would like to emphasize that the matrix-functions 0) and $°"(r], 0), which define 
through the Borel decomposition solutions of 2D Toda lattice, can be written explicitly. 
From (|5.5| ) it follows that <$> a (r),0) is a lower-triangular matrix 

%(v) = 0, i<j, ^(77,0) = 1, (5.15) 

Using ( |5. 15] ) as initial condition for integration of equation ( |5.3|) , c^jy = j, we obtain 

that for i > j entries of $ CT are given by multiple integrals 



*ufa)=/ / •••/ nW^-i)^+fc-i- (5-16) 

** «0 J ^ 



In the same way we may right explicit formulae for the matrix $(£). It is upper triangular 
and if we introduce functions 

q n (0 = exp (-/^n(0^') , &n = e^" +1 

then 

^(0 = 0, i>j, $ii(0 = e- ft(0 , (5.17) 
and for i < j its entries are given by multiple integrals 

= f ...f 3 ff fei-fc^-fc+O^-fc+i. (5-18) 

Let M( n \£,r)) be (n + 1) x (n + 1) matrices which are upper left blocks 

= M ij: 0<ij<n, (5.19) 

of the matrix 

M(e,r ? ) = $(e,0)$ CT (r / ,0)- 1 . (5.20) 
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Then well-known formula for diagonal matrix of the Borel decomposition ( |5.8| ) 

H n = D n = det MW, D. x = 1 

together with Theorem 5 provide an explicit formula for the solution 

normalized by the condition 

?o(0,7?) = 

Note that equations of 2D Toda lattice are invariant with respect to the transformation 
Qn Qn + f{v)i where f(rj) is an arbitrary function. This transformation does not change 
initial data on characteristics. 
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